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Ph, abstract 

Q ' We present a detailed analysis of rcdshift-space distortions in the two-point correla- 

i-i I tion function of the 6dF Galaxy Survey (6dFGS). The i^T-band selected sub-sample 

■ which we employ in this study contains 81971 galaxies distributed over 17000deg^ 

d ' with an effective redshift Zcs = 0.067. By modelling the 2D galaxy correlation func- 

tion, C(rp,7r), we measure the parameter combination f{zcs)as{zes) = 0.423 ± 0.055, 
I ■ where / ~ ^Zii^) ^^e growth rate of cosmic structure and as is the r.m.s. of matter 

^ I fluctuations in 8h~^ Mpc spheres. 

Ijr^ . Alternatively, by assuming standard gravity we can break the degeneracy between 

fSJ ' CTg and th e galaxy bia s para meter, b. Combining our data with the Hubble constant 

; prior from lRiess et all (|2011[) . we measure as = 0.76 ± 0.11 and = 0.250 ± 0.022, 

I consistent with constraints from other galaxy surveys and the Cosmic Microwave Back- 

.,— ^ ' ground data from WMAP7. 

, Combining our measurement of fas with WMAP7 allows us to test the cosmic growth 

^SJ ■ history and the relationship between matter and gravity on cosmic scales by constrain- 

T-H ' ing the growth index of density fluctuations, 7. Using only 6dFGS and WMAP7 data 

ILJ . we find 7 = 0.547±0.088, consistent with the prediction of General Relativity. We note 

. y—{ • that because of the low effective redshift of 6dFGS our measurement of the growth 

I rate is independent of the fiducial cosmological model (Alcock-Paczynski effect). We 

^ . also show that our conclusions are not sensitive to the model adopted for non-linear 

redshift-space distortions. 

Using a Fisher matrix analysis we report predictions for constraints on fas for the 
WALLABY survey and the proposed TAIPAN survey. The WALLABY survey will be 
able to measure fas with a precision of 4 — 10%, depending on the modelling of non- 
linear structure formation. This is comparable to the predicted precision for the best 
redshift bins of the Baryon Oscillation Spectroscopic Survey (BOSS), demonstrating 
that low-redshift surveys have a significant role to play in future tests of dark energy 
and modified gravity. 

Key words: cosmology: observations, cosmological parameters, large-scale structure 
of Universe, surveys, galaxies: statistics 
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The distribution of matter in the cosmos depends on the 
gravitational interaction and the expansion history of the 
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Universe. Assuming that galaxies trace the mass distribu- 
tion, a measurement of galaxy clustering can be used to 
derive fundamental properties of the Universe. 

On large scales the movement of galaxies is dominated 
by the Hubble recession, while on small scales the gravi- 
tational field introduces so-called peculiar velocities. Indi- 
vidual galaxy redshifts combine both Hubble recession and 
peculiar velocities indistinguishably. However, these effects 
can be statistically distinguished in a large sample of galaxy 
redshifts. This is the purpose of this paper. The difference 
between the redshift-inferred distance and the true distance 
is known as redshift-space distortion. Redshift-space distor- 
tions effectively couple the density and velocity fields, com- 
plicating the models needed to accurately describe observed 
galaxy samples. On the other hand they permit measure- 
ments of the properties of the galaxy velocity field, which 
are difficult to access otherwise. In standard gravity we can 
use the amplitude of peculiar velocities to measure param- 
eters that describe the matter content of the Universe such 
as and ug. 

In this paper we report measurements of the parame- 
ter combination /(2;cff)c"8(2cff), where / = dln(Z))/dln(a) is 
the growth rate of cosmic structure (in terms of the linear 
growth factor D and cosmic scale factor a) and erg is the 
r.m.s. of the matter fluctuations in spheres of 8/i~^ Mpc. 
The measurement of /(zefi)crg(2:cff) can be used to test the- 
ories of dark energy and modified gravity, since a stronger 
gravitational interaction causes a larger growth rate /. It 
is interesting to note in this context the fact that a differ- 
ent form of gravitational interaction on large scales could 
be respons ible for the accelerating expansion of the Un i- 
verse (e.g. iDvah. Gabadadze fc Porratil I2OO0I : IWang||2008h . 
Probes such as type la supernovae, the Cosmic Microwave 
Background (CMB) or Baryon Acoustic Oscillations, which 
have proven the existence of th e current accelerati on of the 
expansion of the Universe (e.g. iBlake et al.||201fd '). cannot 
distinguish between acceleration due to a dark energy com- 
ponent with negative pressure or due to a modification of 
General Relativity. However, measurements of the growth 
of structure are able to distinguish between these models. 

In order to measure /(Zeff)o"g(-2cfi) we have to model the 
effect of redshift-space distortions on the correlation func- 
tion. While on large scales linear theory can be used to 
model these effects, on smaller scales non-linear contribu- 
tions complicate the process. Several new approaches have 
been suggested in recent years to extend linear theory. We 
will discuss some of these models and apply them to our 
dataset. 

Redshift-space distortions have previously been anal- 
ysed using both the correlation function and power 
spectrum using da ta from the 2 dF Galaxy Redshift 
Surv ey (2dFGRS ; Feacock et all I2001I : iHawkins et al] 



I2OO3I : ICole et all l2005h and 



Survev (SPSS: iTeginark et al., ..... 
Tegrnark et"ai]|2006l: ICabre fc Gaztanagall2009l : ISong et all 



20fOl : 



Samushia et al 



the Sloan 
2004 



Digita l Sk; 



Zehavi et al. 200 



■ky 



More recently it has become 
possible to do similar stud ies at highe r redshift using the 
VVDS (|Guzzo et all I2OO8I). WiggleZ jiBlake et all l201Ial ') 
and VLT VIMOS surveys (jBielbv et al.lboiol '). 

Redshift-space distortion measurements are also sensi- 
tive to the overall amplitude of the cl ustering pattern o f 
matter, commonly parameterised by ag (|Lahav et al.ll2003 ). 



This parameter is used to normalise the amplitude of clus- 
tering statistics such as the correlation function, ^ oc ag. 
From the CMB we have a very accurate measurement of the 
matter fluctuations in the early universe (the scalar ampli- 
tude As) at the time of decoupling, z*. In order to derive 
(Js{z=Q) we have to extrapolate this measurement to redshift 
zero, involving assumptions about the expansion history of 
the Universe. The CMB constraint on trg heavily depends 
on these assumptions. Hence there is a clear advantage in 
obtaining low-redshift measurements of this parameter. The 
6dF Galax;y Survey, which we analyse in this study, is one of 
the largest galaxy redshift surveys available. Its very small 
effective redshift and wide areal coverage (41% of the sky) 
make it a powerful sample for the study of the local galaxy 
distribution. 

Galaxy surveys usually have to consider degeneracies 
between redshift-space distortions and the Alcock-Paczynski 
effect, which arises from the need to assume a cosmological 
model to transform redshifts into distances. At low redshift 
this effect is very small, meaning that our measurement is 
fairly independent of the choice of the fiducial cosmological 
model. 

While at high redshift {z > f ) the matter density dom- 
inates both the expansion of the Universe and the growth of 
perturbations, at low redshift these two are partially decou- 
pled, with dark energy mostly dominating the background 
expansion and the matter density dominating the growth 
of perturbations. As a result in ACDM and most proposed 
modified gravity models, low redshift measurements of the 
growth rate have a better constraining power than high red- 
shift measurements. The measurement of the growth rate in 
6dFGS therefore not only provides a new independent data 
point at very low redshift, but also promises to make a valu- 
able contribution to tests of General Relativity on cosmic 
scales. 

The outline of this paper is as follows: In section [2] we 
introduce the 6dF Galaxy Survey. In section [3] we describe 
the details of the correlation function estimate and introduce 
the 2D correlation function of 6dFGS. In section |4] we de- 
rive the covariance matrix for the 2D correlation function. In 
section [5] we summarise the theory of redshift-space distor- 
tions, including extensions to the standard linear approach. 
We also discuss wide-angle effects and other systematics, 
such as the Alcock-Paczynski effect. In section |6] we fit the 
2D correlation function to derive gg{zcfi) = /(^cff )o"g(zcff) 
and CTg. Cosmological implications are investigated in sec- 
tion [T] In section |8] we make Fisher matrix predictions for 
two future low redshift galaxy surveys, WALLABY and the 
proposed TAIPAN survey. We conclude in section [51 

Throughout the paper we use r to denote real space 
separations and a to denote separations in redshift-space. 
Our fiducial model assumes a fiat universe with = 0.27, 
w^'^ = -I and 0.1^ = 0. The Hubble constant is set to 
Ho = lOO/ikm s^^Mpc"^ 



2 THE 6dF GALAXY SURVEY 



The 6dF Galaxy Survey (6dFGS: [jones et al.|[20o3 . l2006l . 
l2009h is a near-infrared selected {JHK) redshift survey of 
125 000 galaxies across four-fifths of the southern sky, with 
secondary samples selected in 6j and rp. The |b| < 10° re- 
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Figure 1. The solid black line shows the 6dFGS redshift dis- 
tribution, while the dashed black line shows one of the random 
mock catalogues containing the same number of galaxies. The 
blue solid and dashed lines show the distribution after weighting 
with Po = 1600/i3 Mpc-3 (see section IXTI for more details on the 
employed weighting scheme). 



gion around the Galactic Plane is avoided by the JHK sur- 
veys to minimise Galactic extinction and foreground source 
confusion in the Plane (as is |6| < 20° for bj and rp). 
The near-infrared photometric selection was based on to- 
tal magnitudes from the Two-Micron AH- Sky Survey - Ex - 
tended Source Catalog (2MASS XSC: Ijarrett et al.|[200ol ). 
The spectroscopic redshifts of 6dFGS were obtained with 
the Six-Degree Field (6dF) multi-object spectrograph of the 
UK Schmidt Telescope (UKST) between 2001 and 2006. 
The effective volume of 6dFGS is abo ut the same as the 
2dF Galaxy Redshift Survey (2dFGRS: ICoUess et aLlbOOlh 
and is a little under a third that of the Sloan Digital Sky 
Survey main spec troscopic sample at it s Seventh Data Re- 
lease (SDSS DR7: lAbazaiian et al.ll2009l ). A subset of early- 
type 6dFGS galaxies (approximately 10 000) have measured 
line-widths that will be used to derive Fundamental Plane 
distances and peculiar motions. 

The 6dFGS /iT-selected sample used in this paper con- 
tains 81 971 galaxies selected to a faint limit of 7^ = 12.75. 
The 2MASS magnitudes are on the Vega system. The mean 
completeness of the 6dFGS is 92 percent and median red- 
shift is 2 = 0.05. Completeness corrections are derived by 
normalising completeness-apparent magnitude functions so 
that, when integrated over all magnitudes, they equal the 
measured total completeness on a particular patch of sky. 
This proc edure is out li ned in the luminosity function eval- 
uation of Ijones et al.l (|2006l ) and also in Jones et al., (in 
prep). The original survey papers (jjones et al ] |2004l2009l ) 
describe in full detail the implementation of the survey and 
its associated online database. 

The clustering in a galaxy survey is estimated relative to 
a random (unclustered) distribution which follows the same 
angular and redshift selection function as the galaxy sample 
itself We base our random mock catalogue generation on the 
6dFGS luminosity function, where we use random numbers 
to pick volume-weighted redshifts and luminosity function- 
weighted absolute magnitudes. We then test whether the 



redshift-magnitude combination falls within the 6dFGS K- 
band faint and bright apparent magnitude limits (8.75 ^5 
K s; 12.75). 

Figure [1] shows the redshift distribution of the 6dFGS 
Jf-selected sample (black solid line) compared to a mock 
catalogue with the same number of galaxies (black dashed 
line) . 



3 CORRELATION FUNCTION 
MEASUREMENT 

We calculate the co-moving distances for each galaxy using 
the measured redshift 



Dc = 



Jo 



with 



E{z) = [l 



f^™^(i + z: 



dz' 



1/2 



(1) 



(2) 



where we assume a flat universe with Oj, =0 and fi^ = 
1 — njjf and describe dark energy as a cosmological constant 
(u;""^ = -1). Given the low redshift of our dataset, these 
assumptions have a very small impact on our final results 
(see section [53} . 

We define the positions of two galaxies as si and S2 ■ The 
redshift-space separation is then given by /i = si — S2, while 
s — {si + s*2)/2 is the mean distance to the pair of galaxies. 
Now we can calculate the separation along the line-of-sight 
n and the separation across the line-of-sight Vp 



\s- h\ 



rp=\l\h? 



(3) 
(4) 



The absolute separation is then given by s = y^Tp + rJ. 

We measure the separation between all galaxy pairs in 
our survey and count the number of such pairs in each sepa- 
ration bin. We do this for the 6dFGS data catalogue, a ran- 
dom catalogue with the same selection function, and a com- 
bination of data-random pairs. We call the pair-separation 
distributions obtained from this analysis step DD, RR and 
DR, respectively. In the analysis we used 30 random cat- 
alogues with the same size as the real data catalogue and 
average DR and RR. The redshift-spac e corre lation function 
itself is then given bv lLandv fc Szalavl (|l993l ): 



_ furY ^DR fUr 

where the ratio rir/nd is given by 



n,- 
rid 



(5) 



(6) 



and the sums go over all random (Nr) and data (Nd) galax- 
ies. The galaxies are weighted by the inverse completeness 
d of their area of the sky 



Wi{z) — Ci. 



(7) 



We will discuss further weighting techniques in the next sec- 
tion. 
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3.1 Density weighting 




In Fourier space the error in measuring the amplitude of a 
mode of the linear power spectrurrQ is given by 



Figure 2. The 2D correlation function of 6dFGS using a density 
weighting with Pq = 1600/i'^ Mpc^"^. For reasons of presentation 
we binned the correlation function in 0.5/i~^ Mpc bins, while in 
the analysis we use larger bins of 2h~^ Mpc. Both redshift-space 
distortion effects are visible: the "finger-of-God" effect at small 
angular separation rp, and the anisotropic (non-circular) shape 
of the correlation function at large angular separations. 



There is a possible bias in the estimation of the correla- 
tion function due to the fact that we estimate both the mean 
density and the pair counts from the same survey. This leads 
to a non-zero difference between the true correlation func- 
tion estimate of an ensemble of surveys and the ensemble 
average of ^(s) from each surv ey. This is co mmonly known 
as the integral constr aint (e.g. Peebles! 1980l ). which can be 
calculated as (see e.g. iRoche et al.ll2002 ) 



model RR 

Y.RR 

and enters our correlation function estimate as 

Cdata ~ Cdata + *C, 



(8) 



(9) 



where ^^ata is the redshift-space correlation function from 
eq. [5] and ^modoi is the model for the correlation function. 
In 6dFGS ic is typically around 6 x 10~* and so has no 
significant impact on the final result. 

In Figure [2] we show the 2D correlation function calcu- 
lated from the 6dFGS dataset. In this Figure we use bins 
of Q.hh~^ Mpc, while for the analysis later on we use larger 
bins of 2/i~^ Mpc (see Figure [6|. The figure shows clearly 
the two effects of redshift-space distortions which we will 
discuss later in section [SI the "finger-of-God" effect at small 
rp, and the linear infall effect at larger which gives the 
correlation function a non-circular shape. 



cjp,^k) = {b + fiJ?fP{k) + {N), 



(10) 



where b is the linear bias, / is the growth rate, ^ is the 
cosine of the angle to the line of sight and P{k) is the matter 
power spectrum. The first term on the right hand side of 
this equation represents the sample- variance error, while the 
second term {{N)) represents the Poisson error. 

If the sample-variance error is dominant we can reduce 
the power spectrum error by employing a weighting scheme 
which depends u pon the galaxy density n(z) , such as the 
one suggested bv lFeldman. Kaiser fc Peacockl (| 19941 ) 

"^'(^^= l-fnWo ' 
where Po describes the amplitude of the weighting. A 
stronger weighting (larger value of Po) yields a smaller 
sample-variance error since it increases the survey volume 
by up-weighting sparsely sampled regions. However, such a 
weighting scheme also increases the Poisson error because it 
shifts the effective redshift to larger values with a smaller 



galaxy number density. This is illustrated in Figure 3(a 



and 3(b) Such a weighting scheme is standard for large scale 



structure analyses. 

In a magnitude-limited sample such as 6dFGS, up- 
weighting higher redshift galaxies also has the effect of shift- 
ing the galaxy bias to larger values. The sample- variance 
error is proportional to the clustering amplitude, and so a 
larger bias results in a larger error. However, the weight- 
ing will still ensure that the relative error of the power 
spectrum, (Jp(k)/ P{k), is minimised. The redshift-space dis- 
tortion signal is inversely proportional to the galaxy bias, 
P ~ Q.'^{z)/h. If weighting increases the bias h, it also re- 
duces the signal we are trying to measure. We therefore must 
investigate whether the advantage of the weighting (the re- 
duced relative error) outweighs the disadvantage (increasing 
galaxy bias). 

The situation is very different for measuring a signal 
that is proportional to the clustering amplitude, such as the 
baryon acoustic peak. In this case the error and the sig- 
nal are proportional to the bias, and so weighting will al- 
ways be beneficial. We stress that an increasing bias with 
redshift is expected in almost all galaxy redshift surveys. 
Therefore redshift-space distortion studies should first test 
whether galaxy weighting improves the measurement. The 
6dF Galaxy Survey is quite sensitive to the weighting scheme 
employed because it has a high galaxy density, making the 
sample- variance error by far the dominant source of error. 

Finally, we have to consider the correlation between the 
bins in the measured power spectrum or correlation func- 
tion. If the error is sample- variance dominated, the bins will 
show large correlation (especially in the correlation func- 
tion), while in the case of Poisson-noise dominated errors, 
the correlation is much smaller. Weighting will always in- 
crease the Poisson noise and hence reduce the correlation 
between bins. 

^ As the correlation function and power spectrum are related by 
a Fourier transform, the following discussion also holds true for a 
correlation function measurement. 
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All of the density weighting effects discussed above need 
to be considered, before deciding which weighting is best 
suited to the specific analysis. We summarise as follows: 

• If a galaxy sample is sample-variance limited, density 
weighting will reduce the (relative) error of the clustering 
measurement. 

• In most galaxy redshift surveys, the density weighting 
increases the galaxy bias, which reduces the redshift-space 
distortion signal, by flattening the clustering anisotropy. 

• Galaxy weighting also reduces the correlation between 
bins in the power spectrum and correlation function. 

Whether the signal-to-noise is actually improved by den- 
sity weighting depends on the specific sample. For the 
6dFGS redshift-space distortion analysis we found Po ~ 
1600/i~'^ Mpc^ leads to the most accurate constraint on the 
growth rate. 

This discussion also indicates that low-biased galaxy 
samples have an advantage over a highly biased galaxy sam- 
ple in measuring redshift-space distortions. We will discuss 
this point further in section [S] 



4 ERROR ESTIMATE 

In this section we will derive a covariance matrix for the 2D 
correlation function using jack-knife re-sampling. We also 
use log-normal realisations to test the jack-knife covariance 
matrix. 

4.1 Jack-knife re-sampling 

We divide the dataset into A'' — 480 subsets, selected in R.A. 
and Dec. Each re-sampling step excludes one subset before 
calculating the correlation function. The N — 1 remaining 
subsets have a volume which is (A — 1)/A times the volume 
of the original data. The covariance matrix is then given by 

a, = ^^^^E , (12) 

fc=i 

where £,''{si) is the correlation function estimate at separa- 
tion Si with the exclusion of subset k. The mean value is 
defined as 

?(^') = ^E«'(^0- (13) 

The case i = j gives the error ignoring correlations between 
bins (Ti = Cii. 

4.2 Log-normal realisations 

We can create a l o g- normal r ealisation jColes fc Joneslll991 ; 
Cole et all l2005l: jpercival et al., .20071 : iBlake et all 120113 : 
Beutler et al.ll201lh of a galaxv survev by deriving a density 
field from a model power spectrum, P{k), assuming Gaus- 
sian fluctuations. This density fleld is then Poisson sam- 
pled, taking into account the window function and the total 
number of galaxies. The assumption that the input power 
spectrum has Gaussian fluctuations can only be used if the 
fluctuations are much smaller than the mean density, oth- 
erwise the Gaussian model assigns a non-zero probability 




7i[h-'Mpc] (r^ = 11h-'Mpc) 

(a) 



10 F 

- Poisson: = 20 000 h'^Mpc^ 

o total: P, = 20 000 h 'Mpc' 

2 Poisson: = 1 600 ti Mpc^ 

• totai: Pq = 1 600 li'^iVlpc^ 




5 10 15 20 25 30 

Tp [h-^Mpc] (ii = 11IT'Mpc) 

(b) 

Figure 3. (a) The relative error in the 2D correlation function as 
a function of line-of-sight separation tt at a fixed rp = llh~^ Mpc. 
Other regions of the 2D correlation function behave in a similar 
manner. The solid lines show the Poisson error for different values 
of Pq, while the data points show the total (Poisson -|- sample vari- 
ance) error obtained as the diagonal of the covariance matrix de- 
rived using jack-knife re-sampling. The purpose of the weighting 
(Po) is to minimise the total error, which is achieved for a value of 
Po ~ 1600/i~'' Mpc'^. The weighting reduces the error by almost 
a factor of two on most scales. The dashed line shows the error 
derived from log-normal realisations using Po = 1600h-3 Mpe^ 
and is in very good agreement with the jack-knife error, (b) same 
as (a) for a fixed line-of-sight separation tt = llh~^ Mpc. 

to regions of negative density. A log-normal random fleld, 
LN{x), can avoid this unphysical behaviour. It is obtained 
from a Gaussian field G{x) by 

LN[x) = exp[G(f)] - 1, (14) 

which is positive-definite but approaches G(x) whenever the 
perturbations are small. As an input po wer spectrum for the 
log-normal field we use the linear model (|Kaiserlll987h . given 
by 

P,{k,^JC)^h\l + P^^Pssik) (15) 

with b — 1.47 and /? = 0.35, which is consistent with the 
parameters we measure for 6dFGS (see section (6)1 . Pss (fc) 
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Figure 4. The correlation matrix for the 2D correlation function 
^(rp,7r) with a bin size of 2 X 2/i^^Mpc. The upper-left corner 
shows the jack-knife estimate, while the lower-right corner shows 
the result of using 1500 log-normal realisations. Since this plot 
shows the correlation of all 15 X 15 bins it contains 225 X 225 
entries. 
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Figure 5. This plot shows the correlation of bin 127 (rp = 
Vihr^ Mpc, TT = nh^^ Mpc) with all other bins in the 2D correla- 
tion function, derived using jack-knife re-sampling. It corresponds 
to row/column 127 of the jack-knife correlation matrix which is 
shown in Figure |4] (upper- left corner). 



is a linear de nsity power spect rum in real space obtained 
from CAMB (|Lewis et al.ll2000l ') and Pg{k,fi) is the galaxy 
power spectrum in redshift-space. Ou r method is explained 
in more detail in lBeutler et al.l (|201lh . appendix A. 

We produce A'' — 1500 such realisations and calculate 
the 2D correlation function for each of them, deriving a co- 
variance matrix 



1 ^ 

Cij = Jj—TJ [^"(f"?. 1") - ^'{rp, 7r)j 

n = l 

X [Cn(rp,7r)-C^(rp,7r)] , 
where the mean value ^^('"pi'''") is defined as 
— 1 ^ - 



(16) 



in rough agreement. The diagonal errors for both the log- 
normal and jack-knife covariance matrices are plotted in Fig- 
ures 3(a) and 3(b) and show very good agreement at large 



scales. 

Every row/column in Figure [3] shows the correlation 
of one bin with all other bins in (,{rp,n). Figure [5] shows 
an example of such a row/column obtained from jack-knife 
re-sampling as a 15 x 15 matrix, in this case for bin 127 
(rp = 13h-^ Mpc, TT = 17/i"^ Mpc). 

Log-normal realisations do not account for non-linear 
mode coupling and are very model-dependent in the quasi- 
linear and non-linear regime. Since our analysis relies on fits 
to fairly small scales, we decided to use the jack-knife co- 
variance matrix in our analysis. However, we find that none 
of the results reported in this paper depend significantly on 
which of the two covariance matrices is used. 



(17) 



and QiirpjTr) is the 2D correlation function estimate of re- 
alisation n at a specific separation (rp,7r). 



4.3 Discussion: Error analysis 

Figure 2] shows the correlation matrix 



(18) 



derived from the two different covariance matrices, dj. We 
plot the jack-knife result in the upper-left corner and the 
log-normal result in the lower-right corner. Both the cor- 
relation and covariance matrix are symmetric. While the 
log-normal correlation matrix indicates somewhat more cor- 
relation between bins, overall the correlation matrices are 



5 MODELLING THE 2D CORRELATION 
FUNCTION 

In this section we discuss the theory of redshift-space distor- 
tions, starting with the standard linear perturbation theory. 
We then describe different approaches for extending the lin- 
ear model to include non-linear structure formation. We also 
discuss deviations from the plane-parallel approximation. 



5.1 Linear redshift-space distortions 

The position of a galaxy in real space r — {x, y, z) is mapped 
to the position in redshift-space, s, via 



s = r H r:^T^ z 



(19) 



aH{z) 
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where the unit vector z indicates the line-of-sight direction, 
and the quantity v- is the line-of-sight component of the ve- 
locity field, namely Vz = v-z. The scale factor a is defined as 
1/(1-1- z) and H{z) is the Hubble constant at redshift z. The 
second term in the equation above represents peculiar ve- 
locities caused by gravitational interaction. On small scales 
this elongates structures along the lin e-of-sight and leads 
to the so-called "finger-of-God" effect (|jackson|[l973) . On 
large scales matter falls in towards over-dense regions and 
systematically influences our distance measurement, making 
the over-densities appear more over-dense. The latter effect 
can be described by linear theory, while the "finger-of-God" 
effect is a non-linear phenomenon. 

The mod el of linear re dshift -space distortion s has been 
developed by iKaiseJ (|l987l ) and iHamiltonl (|l992l ) assummg 
a plane-parallel approximation. In Fourier space the linear 
model can be written as 



(20) 



where Pss{k) is the matter density power spectrum and 
Pg{k, fi) is the galaxy density power spectrum. In this model, 
linear redshift-space distortions are quantified by the param- 
eter P, which is defined as 



b dln(a) b ' 



(21) 



where b is the linear galaxy bias factor, D{z) is the growth 
factor and 7 is the gravitationa l growth ind ex, which takes 
the value 7 = 0.55 in ACDM (|Lindeij|2005l 'l. n^iz) is the 
matter density at redshift z, and is defined as 



nm{z = o){i + zf 



H{Z) 



with 



H{z) 



(22) 



(23) 



where we again follow our fiducial model of £7^ =0 and 
w^"^ = —1. Modifications of the gravitational force mainly 
affect 7, while changes in the expansion history of the Uni- 
verse affect Q.m{z). 



5.3 Extensions to linear theory 

For the remainder of this section we will discuss possible 
extensions of the linear model to include non-linear struc- 
ture formation and wide-angle effects. The simples t model 
of non-linearities is the so- called streaming model (|Peeblesl 
Il98d : iHatton fc Colelll998h . Here, the redshift-space corre- 
lation function is just a convolution of the linear correlation 
function in redshift-space with a pairwise velocity probabil- 
ity density function, F{v), 



Cst(rp,7r) 



H{Ze{{)aeB 



F{v)dv, (26) 



where Oefr = 1/(1-1- Zeff) is the scale factor and H{zcff) is the 
Hubble constant at the effectiv e redshift. We chose F(v ) to 
be an exponential distribution (jPeacock fc Dodd3ll996h 



Fiv) 



exp 



-V2\v\ 



(27) 



which has been shown to succ es sfully descri b e ob- 
servations JPavis fc Pee bles' 'l982|; iFisher et all Il994l : 
iMarzke et"ar.l 1 19951 : iLandv, i2002i ). Although the pairwise 



velocity dispersion within a halo is expected to follow a 
Gaussian distribution instead of an exponential, galaxies 
populate halos of a wide range of masses and velocity 
dispersions, which c ombine to approximately form an 
exponential functi on j|Sheth|[l99^ : iDiaferio fc Gelled Il996l : 
ISeto fc Yokovama 1 19981'). The paramet er ap depends on 
galaxy type (e.g. Maderwick et al.ll2003l ) and hence its use 
for cosmological constraints is limited. 

In recent years many improvements to the model dis- 
cussed above have been suggested. We will initially discuss 
these models in Fourier space because the theoretical moti- 
vation is clearer and many expressions needed for these mod- 
els simplify considerably. However, since our data is in real 
space we will also give real-space expressions later on. We 
start with the streaming model of eq. 1261 which in Fourier 
space is given by 



Pg(k,^)=b^(l + P^^yPss(k) 



(28) 



5.2 Parameterisation 

From now on we will formulate our equations in terms of 
ge{z) — f{z)as{z) and gb{z) — bas{z). We choose the pa- 
rameter set [go, gb] instead of [/3, as, b], because erg and 
the lin ear bias, b, are degenerate and difficult to disen- 
tangle (IWhite. Song fc Percivalll2008l : ISong fc Percivalll2009l : 
ISong et al.ll2010D . Within this parameterisation the power 
spectrum is expressed as 

Pg{k) = b^Pssik) = gtQssik), (24) 

where Qss{k) is the unnormalised matter density power 
spectrum (see eg. I30p . In terms of our new parameters, we 
can write 

/3=^. (25) 
9b 



Hence in Fourier space, the convolution with an exponential 
function becomes a multiplication by a Lorentzian distribu- 
tion. 

The model above assumes that there is a perfect corre- 
lation between the velocity field and the density field, which 
is given by Pss{k) — Pse{k) — Pee{k) where Pss{k) is the 
matter density power spectrum as before. Pge{k) — {\dk\^) 
is the velocity divergence power spectrum where 9 — V -v is 
the velocity divergence, and Pse{k) is the cross power spec- 
trum. Non-linear effects will violate these assumptions, since 
the density power spectrum is expected to increase in am- 
plitude at small scales because of non-linear effects, while 
the velocity field becomes randomised at small scales (e.g. 
within virialized galaxy clusters) and h ence Pee(k) will de - 
crease in amplitude (e.g. ICarlsonll2009l ) . IScoccimarrol (|2004l ) 
suggested expressing the 2D power spectrum without the 
assumption of linear relations between the density field and 
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velocity field, by 

Pg{k,fl) = Fq{k,fj,,a^) 

X [b^Pssik) + 2fi^bfPse{k) + fi\fPee{k)] 
= Fq{k,fJ.,crv) 

X [9bQss{k) + 2^^gi,geQse{k) + ^*glQee{k)\ 
where the different Q^y are defined as 

Qss{k) = P44(fc)/(J8(2cff)^, 

Qse{k) = P5e{k)/as{zcfi)^ , 
Qee{k) = Pee{k)/(Js{zcfi)^ . 



(29) 



(30) 



and the damping function Fq{k, fi, a^) is usually chosen to 
be a Gaussian of the form 



Fq(k, ^, (T„) = e 



(31) 



The parameter (t„ quantifies the non-linear dispersion in the 
bulk motion of halos. It is different to the ap parameter 
we introduced earlier that describes small-scale randomised 
motion (e.g. of galaxies within a halo). We can derive 
from the velocity power spectrum as 



2/ X gejzf 

^"(^^ ^ 2 



9{k)dk. 



(32) 



Ijennings. Baugh fc Pascolil Honi) provide fitting formu- 
lae for Pseik) and Poe(k) derived from N-body simulations. 
They find the following relation between the different power 
spectra 



Pxyik) 



y^P^ + aiPi,{k) 
a2 + a3Pssik) 



(33) 



where Pss(k) can be obt ained from CAMB by including 
halofit (jSmith et all 120031 ;). For the cross power spectrum 
Pxyik) — Pse{k) we use the (updated) parameters {ao, 
ai, 02, Q3) = (-12 483.8, 2.55430, 1381.29, 2.54069) and 
for Pa:y{k) = Peg(k) we use (-12480.5, 1.52404, 2165.87, 
1.79640). The fitting formula reproduces Pee to better than 
1% for k < 0.4/1 Mpc-\ to 10% for 0.4 < fc < 0.7/iMpc-\ 
and to 15% for 0.7 < k < l/iMpc~^. It also reproduces Pse 
to less than 4% over the whole range k < l/iMpc~^ (Jen- 
nings, private communication). We cut off the integral when 
the Jennings formula predicts negative values, although be- 
cause of the high precision of the fitting formula up to large 
k, such a cut-off will not affect our measurement. 
We can express eq. [29] in real-space as 



^Sc(rp,7r) = 



2 1 

lUo.ssir) + i:gbge(.o,se{r) + -glio,ee{r) 



+ 



4 4 2 

-^gbge^2,se{r) + -ge^2,ee{r) 



5- 
^2(m) 



Voifi) 



+ -;^geU,ee{r)'PA(p), 
35 



(34) 



where VeifJ,) are the Legendre polynomials and the spherical 
harmonic moments ^e.xy{r) are given by 



(35) 



X coii{kr^)Qxy(k)'Pi{^). 



Appendix|X]shows a partial analytic solution for the double 
integral above. 

We therefore have two different models which we will 
apply to our data: The simple streaming model £,st{rp,n) 
and the Scoccimarro models 5sc(''p,7r). Note that (,Sc{rp,n) 
does not include the parameter (jp and hence has one less 
free parameter than the streaming model. 

All equations above are based on the plane-parallel ap- 
proximation. In the case of 6dFGS we also need to account 
for wide-angle effects. This means we will replace the equa- 
tions above with more general descriptions, which we discuss 
in the next section. 



5.4 Wide-angle formalism 

So far we have assumed that the separation between galaxy 
pairs is much smaller than the distance of the galaxies 
from the observer. The 6dF Galaxy Survey has a maxi- 
mum opening angle of 180° and the (effective) redshift is 
fairly low at Zeff = 0.067 (see section (3)1. We therefore in- 
clude wide-angle correction terms. The wide-angle descrip- 
tion of r edshift-s pacc distortions has been laid out in several 
papers dSzalav c t ah, , 1997 ; Szap udi 200 4;_Ma tsubara, 20o3 : 
iPapai fc Szapudill2008inRaccaneUi et al.ll2010l 'l. Here we wiU 
expand on this work by formulating the equations in terms of 
ge and gt and by distinguishing between the density-density, 
velocity- velocity and density- velocity contributions {(,ss,(,ee 
and (,se)- The model will then correspond to the Scoccimarro 
model (eq. I34|l in the last section. 

The general redshift-space correlation function (ignor- 
ing the plane-parallel approximation) depends on (j), 6 and 
s. Here, s is the separation between the galaxy pair, 9 is the 
half opening angle, a nd (j) is the angle of s to the line-of- 
sight (see Figure 1 in iRaccanelli et al]|2010l ). The angles (j) 
and 6 are not independent, but the relation between them 
is usually expressed through the two angles 0i and 02 given 
hy (f) = \{4>i + 02) and 6 = \{4>i — 4>2)- The total corre- 
lation function model, including 0(6^ ) correction terms, is 



then given by (|Papai fc Szapudill2008h 



5(0, 6, s) = aoo + 2ao2 cos(20) -I- 022 cos(20) + 622 sin^(20) 
+ [ - 4ao2 cos (20) - 4a22 - 4^22 - 4aio cot^(0) 
-I- 4aii cot^(0) - 4ai2 cot^(0) cos(20) + Abu 



Bl2 COS 



(36) 



This equation reduces to the plane-parallel approximation if 
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0. The factors a^y and bxy in this equation are given by 



5.5 Systematics and the Alcock-Paczynski effect 



aoo = 9b£.o,ss{r) H —£.o,seir) + -j^5o,es('") 



dbge <;2 / N , Sge 2 / x , Sge 2 / s 
-3-42,ie(f-) + —^2.eB{r) + Tj7^U,ee[r) 



21 



140 



ao2 = 5-6, + -T-rG.eel'^) + ir^U.eeir) 



2 

9s f 2 



14 



28 



C122 = — Co.eel'^) - 7^^2,oei^) + -^77^4, eel'') 



21 



140 



h ys t2 / s ye c2 / s it/e ,2 / ^ 
t>22 = —^o.oe[r) - -^^2,oe{r) - -^U.oe{r) 



aio = 



5 r or 



(37) 



ai2 = ^[2eiMr)-^lMr)] 



&12 = 



where the spherical harmonic moments, ^^"^.^(t), are 
X cos{krfi)Qa:y{k)Ve{fJ.) 



(38) 



with Qss{k), Qse{k) and Qee(k) as defined in equation [30] 
(see appendix 1X1 for an analytic solution). 

In order to obtain a model for the 2D correlation func- 
tion (including wide-angle effects), we can simply integrate 
eq. [36] over 9. In case of the plane-parallel approximation 
this would correspond to eq. 1341 To reduce the equation 
to the simple streaming model we have to set cr^ = 0, 
Qss ik) = Qsejk) = Qee(k) and convolve with F{v). 

ISamushia et al.l (|2011 ) studied wide angle effects in the 
SDSS-LRG sample and found that such effects are very 
small. The lower redshift and larger sky coverage of the 6dF 
Galaxy Survey mean that wide-angle effects are certainly 
larger in our data set than SDSS-LRG. However, we found 
no significant impact of such effects on our results. This 
is mainly because our analysis includes only rather small 
scales (tt < 30/i~^ Mpc, Tp < 30fe~^ Mpc, see sect ion [6]). 
This agrees with the findings of iBeutler et al.l ()201ll ). 

We stress that the correction terms discussed above only 
capture first-order effects and so we additionally restrict our 
analysis to ^ < 50°. 

Many authors have found, using N-body simulations, 
that both the simple streaming and Scoccimarro mod- 
els do not capture all non-lin ear and quasi-linear effects 
presen t in redshift surveys (e.g . [Jennings. Baugh fc Pascoli 



[2011bl : [Kwan. Lewis fc LindeJ [2OIII : [Torre fc Guzzd [2oTs 
and references therein). For example the linear bias model 
seems to be too simplistic at small scales. In this study 
we restrict ourselves to the models discussed above but 
would like to point out that m any other models have been 
suggested ([Matsubara[[2008a[lbl: [Taruva. Nishiniichi fc Saitol 



[2OIOI : [Reid fc White|[201l[ : [Seliak fc McDonald|[201ll ) which 
could be compared to our data in future work. 



A dataset such as 6dFGS, containing galaxies with a high 
linear bias factor, may be prone to scale-dependent galaxy 
bias on small scales. We used the GiggleZ simulations (Poole 
et al., in preparation) to derive the form of the scale- 
dependent bias. We find a 1.7% correction to a linear bias at 
scales of s = lOh^^ Mpc, and note that this correction has 
a small but non-negligible impact on our results, depending 
on the smallest scales included in the fit. 

A further systematic uncertainty comes from the so- 
called Alcock-Paczynski (AP) effect ([Alcock fc Paczvnskil 
[l979t ). By assuming a cosmological model for the con- 
version of redshifts and angles to distances (see sec- 
tion [3]) we could introduce an additional anisotropic sig- 
nal in the correlation function which depends on the 
angle to the line-of-sight, if the true cosmology differs 
from our fiducial cosmological model. These distortions 
are degenerate with the linear redshift-space distortion 
signal and hence bot h effects need to be modelled to 
avoid systematic bia s ( Ballinger. Peacock fc Heavens! [l 9961 : 
[Simpson fc Peacock|[2010l ). Using the two scaling factors 



/II 
f± 



H[z) 
Da{z) 



(39) 
(40) 



we can account for the AP effect by rescaling the correspond- 



mg axis 



TT 

h 

f± 



(41) 
(42) 



While this is a matter of concern for high-redshift surveys, 
the 6dF Galaxy Survey is very nearly independent of the 
Alcock-Paczynski distortions because the distances tt and 
Tp axe almost independent of the fiducial cosmological model 
when expressed in units of ft^^Mpc. For example the value 
of /x at z = 0.067 changes by only (0.3%, 1.5%) for 10% 
changes in {flm, w). The corresponding values for /y are 
(0.7%, 1.6%). The relative tangential/ radial distortion de- 
pends on the value of (1 + z)Da{z)H{z)/c relative to the 
fiducial model: 



Da{ z)H{z) 



fj_ 

h' 



(43) 



This parameter combination changes by only 0.3 (0.1)% for 
10% changes in £7™ (ui). 

Finally, it has been shown that 7 has a degener- 
acy wath__tlM_e2uatto parameter for dark energy, 
w ([Simpson fc Peacock[[20l"ol ). High redshift measurements 
usually assume ui = — 1, which could introduce a bias in the 
measured value of 7 if dark energy is not exactly a cosmo- 
logical constant. Again 6dFGS is robust against such effects. 



6 FITTING THE 2D CORRELATION 
FUNCTION 

We now fit the two models for the 2D correlation function 
we developed earlier (5st('"p, ti") and ^Sc('"p, ti")) to our data. 
For the final constraints on /erg we use the £,Bc(rp,Ti) model 
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Table 1. Cosmological parameters derived from the 6dFGS 2D 
correlation function. The effective redsfiift is Zeff = 0.067. The 
last column indicates the priors/assumptions which go into each 
individual parameter measurement. T he prior on the Hubble 
const ant comes from iRiess et al. and the WMAP7 prior 

from iKomatsu et al The asterisks denote parameters 

which are derived from fitting parameters. 



Summary of parameter 


measurements from 6dFGS 




0.423 ± 0.055 




9()(^cff) 


1.134 ±0.073 




/3* 


0.373 ± 0.054 




T8 


0.76 ±0.11 


[Ho = 73.8 ±2.4, 7 = 0.55] 




0.250 ± 0.022 


[Ho = 73.8 ± 2.4,7 = 0.55] 


b 


1.48 ±0.27 


[Ho = 73.8 ± 2.4,7 = 0.55] 


/*(2off) 


0.58 ±0.11 


[ge + (T8 from 6dFGS] 


7 


0.547 ±0.088 


[WMAP7] 




0.271 ±0.027 


[WMAP7] 



since it gives similar results to the streaming model with one 

less free parameter. 

Other studies fe.g. lSamushia et al.ll201ll ) prefer to anal- 
yse the correlation function moments ^o, ^2 (and if possible 
^4), which carry the same information as the 2D correlation 
function. The correlation function moments have the advan- 
tage that the number of bins grows linearly with the highest 
scales analysed, while for £,{rp,Ti) the number of bins grows 
quadrati cally. This makes it eas ier to get reliable covariance 
matrices. ISamushia et al] (|201ll ) also show that the measure- 
ment errors of are more Gaussian. However, the correla- 
tion function moments are integrals over ^ and hence carry 
information from all directions, including ^ = Q. Finger-of- 
God distortions can influence the correlation function mo- 
ments up to large scales (20 — 30/i~^ Mpc), while in the 2D 
correlation function they can be excluded via a cut in r^. 
In 6dFGS we have found that these non-linear effects have 
a strong impact on the correlation function moments up to 
30/i~^ Mpc. We have therefore decided to focus on the 2D 
correlation function instead of the correlation function mo- 
ments. 

6.1 Derivation of the growth rate, gg — fas 

In Figure [6] we show the 6dFGS 2D correlation function. 
For our analysis we bin the data in 2 x 2h^^ Mpc bins from 
to 30/i~^ Mpc in Vp and tt. Including larger scales does 
not add further information. At small Vp, the finger-of-God 
effect becomes dominant and we expect any linear model to 
fail. Since our description of non-linearities, in both of our 
models, is limited in its capability to capture all non-linear 
effects, it is necessary to include a cut-off scale Tp"* marking 
a lower limit of the fitting range in Vp. 

Figure [7] shows the measured value of gg as a function 
of the cut-off scale rp"' for our two different models. Above 
rp"' ~ Sh^^ Mpc the streaming model, ^st(''p, t), approaches 
a constant value of gg. Our second model, ^Sc(rp,7r), con- 
tains a systematic error up to much larger scales, before it 
comes into agreement with the streaming model at about 
Tp"' — I6h~^ Mpc. This is expected since this model does 
not include a description of effects in the non-linear regime. 




Figure 6. The 2D correlation function in 2h~^ Mpc bins. The 
fitting area is indicated by the dashed lines, where black corre- 
sponds to the streaming model, (Cst(''p, tt)) and red corresponds 
to the Scoccimarro model (Csc('"p! tt)). The black and red contours 
show the best fitting models for ^st(''p, f) (black) and Csc(''pi''") 
(red). The deviations seen in the two contours at large scales 
are well within the error bars of the two models, which can be 
seen in Figure[8] At small scales (< 14/i~^ Mpc) the Scoccimarro 
model predicts much more clustering, while in the real data this 
clustering is smeared out along the line of sight because of the 
finger-of-God effect. 



For the final constraints we choose Tp"' = lOh^^ Mpc for the 
streaming model and rp"* = 16/i~^ Mpc for the Scoccimarro 
model. We also note that since the Scoccimarro model is 
based on only two free parameters {ge and gb), the error is 
generally smaller compared to the streaming model, which 
has three free parameters (ge, gt and CTp). Oth er studies fit 
for the parameter cr^ (e.g. iTorre fc Guzzd[2012l ) in the Scoc- 
cimarro model, but we derive it using eq. 1321 

For ^sc{rp,Tv) we use the fitting range < tt < 
30h~^ Mpc and 16 < rp < 30h~^ Mpc, which results in a to- 
tal of 105 bins. The best-fitting results are gg = 0.423±0.055 
and gt — 1.134 ± 0.073, where the errors for each parame- 
ter are derived by marginalising over all other parameters. 
The of this fit is 115 with 103 de grees of freedom (d.o.f.), 
indicating a good fit to the data. 

For (,st{rp,TT), we have the fitting range < tt < 
30/i"^ Mpc and 10 < rp < 30fe"^ Mpc, which results in 
a total of 150 bins. The best fitting parameters are gg — 
0.389±0.067, gt = 1.084±0.036 and o-p = 198±81 km/s. The 
reduced of this fit is given by /d.o.f. = 202/147 = 1.37. 
We compare the constraints on gg and gb from both models 
in Figure [8] 

In the Scoccimarro model we could use the param- 
eter (j„ PC gg instead of gg to test cosmology, as sug- 
gested bv lSong et"aLl (|2010t ). Our best fit gave = 2.59 ± 
0.34/i~^ Mpc. However, this parameter depends on an addi- 
tional integral over the velocity power spectrum, which adds 
a theoretical uncertainty. We therefore prefer to use gg in the 
following discussions. 
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Figure 7. The value of ge(zoff ) = /(^cff )''"8{-2cff ) a-s a function of 
the cut-off scale r^j"', obtained by fitting the 6dFGS 2D correla- 
tion function with two different models (as described in section fS.SI 
and l5.4|l . At large scales the two models converge to similar val- 
ues, while on small scales the models deviate from each other 
because of the different descriptions of non-linear evolution. For 
the final parameter measurements in Table [l] we chose model 2, 
?Sc('"pi "■), with a conservative cut-off scale of r^"' = 16/i^^ Mpc. 
fn the lower panel we plot the reduced as an indicator of the 
quality of the fit. 



Figure 8. Likelihood distribution of gg and gi, derived from the 
fit to the 2D correlation function. The solid black contours show 
model ^Sc (''p , I") , while the dashed contours show the streaming 
model (see section 15.31 and 15.41 for details of the modelling) . The 
fitting range is < tt < 30/i~^ Mpc and 10 < rp < 30/i~^ Mpc 
for ^stC^p, tt) and < tt < 30h"^ Mpc and 16 < r-p < 30/1"^ Mpc 
for CscC^pi"")- The black cross indicates the best-fitting value for 
the solid black contours. 



We can also express our results in terms of j3 which 
is given by /? = ge/gb ~ 0.373 ± 0.054. We summarise all 
measured and derived parameters in Table [1] 

6.2 Derivation of erg and fim 

In this section we use redshift-space distortions to directly 
measure ag. The angular dependence of the redshift-space 
distortion signal in the 2D correlation function allows us 
to measure /3, which quantifies the amplitude of redshift- 
space distortions. Together with Qrn.{z) and 7 = 0.55, this 
constrains the linear bias b through the equation 

6~^^^. (44) 

Knowing 6 we can use the absolute amplitude of the 
correlation function, [bas{z)]'^, to constrain ag^z—Q) = 
[D{z=0)/D{z,a)] X asM. 

For computational reasons we use our first model, 
^st {rp, it), in this sub-section and fit the five parameters 
erg, Qrn, b, Ho and ap using an MCMC approach. Since 
the shape of the correlation function is only sensitive to 
r = Qmh, we cannot constrain Q,rn. and Ho at the same 
time. For the final results we include a pr ior on the Hubble 
constant {Ho = 73.8 ± 2.4kms~^ Mpc~\ lRiess et al.|[2oTll . 
from now on referred to as HST prior) and marginalise over 
it. We use the same binning and fitting ranges as in the 
previous section. 

The best-fitting model results in x^/d.o.f — 1.35. We 
find erg = 0.76 ± 0.11, = 0.250 ± 0.022, b = 1.48 ± 0.27 




I I I I \ I I I \ I I I I I I I L_ 
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Figure 9. This plot shows the likelihood distribution of the 
galaxy bias b and eg, which we obtained by fitting the 6dFGS 
2D correlation function assuming 7 = 0.55. The solid black 
line shows the result using a prior on the Hubble c onstant of 
Ho = 73.8±2.4kms-l Mpc-i from lRiess et al.l boill) . while the 
dash ed black line uses a prior of Ho = 67 ± 3.2kms~^ Mpc~^ 
from iBeutler et al.l 1I2OIII ). Although the detection of redshift- 
space distortions can partially break the degeneracy between b 
and erg which exists in the ID correlation function, there is still a 
significant residual degeneracy. The black cross marks the maxi- 
mum likelihood value for the solid black lines. 
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and Op = 174±73km/s. The remaining degeneracy between 
the bias h and ag is illustrated in Figure (9] We include all 
these results in Table [T] 

Figure [TUl compares the 6dFGS f2m — crs probability dis- 
tribution to measurements from several other data sets: The 
CFH T wide synoptic Legacy Survey (CF HTLS) (IFu et al 



200S), the SFI-I— |- peculiar velocity survey jNusser fc Davi: 
2011), cluster abund ance from X-ray surv eys (jMantz et al 
20101 ) and WMAP7 (|Komatsu et alil201lh . Many of the ex- 



periments shown in this figure have systematic modelling 
uncertainties when extrapolating to 2 = arising from 
assumptions about the expansion history of the Universe. 
Only 6dFGS and SFlH — h are at sufficiently low redshift to 
be independent of such effects. To illustrate the impact of 
these effects we plot in Figure[TT]the probability distribution 
f2m — o"8 from WMAP7 for different cosmological models. 
The CMB measures the scalar amplitude As, which needs 
to be extrapolated from redshift z* ~ 1100 to redshift zero 
to obtain erg . Every parameter that influences the expansion 
history of the Universe in this period affects the value of erg 
derived from the C MB alone. The relation between As and 
(Tg is given by (e.g. iTakada. Komatsu fc Futamasell2004l ) 



As 



2c' 



3sin(fc7?) 3cos(fci?) 
{kRy 



(fci?)2 



(45) 



where R = 8 fe~^ Mpc, k, = 02M pc~^ and As = (2.21 ± 
0.09) X 10"^ (jKomatsu et al.ll2009h . D{k,z) is the growth 
factor at redshift z and T{k) is the transfer function. If 
higher-order cosmological parameters such as the dark en- 
ergy equation of state parameter w are marginalised over, 
then the measurements of Qrn and as weaken considerably 
(see Figure [TT|| . 

We now asses s the influence of the Ho prior. We re- 
place the result of iRiess et al.l |2oi3) with a measurement 
derived fr om the 6dFGS dataset using Baryon Acoustic Os- 
cillations (jBeutler et al.|[201ll ). The prior from this study 
is lower than the former value and is given by Hq = 
67±3.2 kms~^ Mpc"^ Using the 6dFGS value of Ho results 
in erg = 0.75±0.13, fi™. = 0.279±0.028, b = 1.52±0.29 and 
o-p = 174 ± 106km/s. The quality of the fit is x^/d.o.f. = 
1.35, very similar to the value obtained with the HST prior. 
Comparing the two results shows that a different prior in Ho 
shifts the constraint in erg and b along the degeneracy shown 
in Figure |9] However, we note that the 6dFGS measurement 
of Ho is derived from the same dataset as our present study 
and hence could be correlated with our measured growth 
rate. We use these results only for comparison, and include 
the values obtained using the HST prior in Table [T] as our 
final results of this section. 

Alternative methods for deriving erg or fag at 
low redshift are provided by peculiar velocity surveys 
fe.g. iGordon. Land fc Slosaj l2007l: lAbate fc Erdogul |2009|; 



Nusser fc David I2OIII: iTurnbuU et al.1 I2OIII : iDavis et al] 
2OIII : [Hudson fc Turnbullll2012l ). 6dFGS will soon provide 
its own peculiar velocity survey of around 10 000 galaxies. 
Velocity surveys have the advantage of tracing the matter 
density field directly, without the complication of a galaxy 
bias. However, they are much harder to obtain and current 




CFHTLS 
(L. Fu 2008) 
SFI++ 

(A. Nusser 2011) 
X-ray cluster 
(A. Mantz2010) 
WMAP7 

(E. Komatsu 2010) 



Figure 10. This plot shows the likelihood distribution in erg 
and Qrn for different cosmological probes. The solid blue con- 
tours show the 6dFGS result, the magenta dotted dashed con- 
tours shows the currently best result of weak lensing from the 
CFHT wide synoptic Legacy Survey i IFu et aLlliool) . the red 
soUd contours show the res ult of the SFlH — h peculiar velocity sur- 
vey llNusser fc D avis'201ll). the green dashed contours show the 
result of iMantz ct al. (201o|) using cluster abundances and the 
black dotted contours arc from WMAP7 jKomatsu et al.ll201ll ). 
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Figure 11. Likelihood distribution of erg — Hm from WMAP7 
for a ACDM model (solid black line), «)CDM model (dashed red 
line), oCDM model (dotted dashed green line) and ouiCDM model 
(dotted magenta line). In blue we show the 6dFGS result. Both 
parameters are defined at redshift zero and WMAP constraints 
on these parameters depend on assumptio ns about the expansio n 
history of the Universe. We use CosmoMC llLewis fc Bridlel2002^ , 
together with the WMAP7.2 l|Komatsu et al.ll201ll ') dataset to 
produce these likelihood distributions. 
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Figure 12. Comparison of measurements of the growth of struc- 
ture using galaxy surveys at different redshifts. The different 
data points belong to 6dF GS (solid blue circle , this paper), 
2dFGRS (so lid blaclc circle; [Hawkins et al.ll2003h . SDSS (solid 
blac k boxes: [ Sa rnushia et alj 201ll ). WigeleZ (solid black trian- 
gles: iBlake etal.l 120111 7 and VVDS (empty circle; IGuzzo et al.l 
l2008t) . We also included a WALLABY forecast with a 4% error- 
bar in red (sec secti on [Si. For the Wiggl eZ survey we also include 
the data points from lBlake et al] ll2011lj ) (empty triangles, shifted 
by Az = 0.005 to the right for visibility), where the Alcock- 
Paczynski effect has been taken into account. We plot a ACDM 
model as well as a DGP model for comparison. 



velocity surveys are 1 — 2 orders of magnitude smaller than 
galaxy redshift surveys. 



7 COSMOLOGICAL IMPLICATIONS 

In this section we test General Relativity by mea- 
suring the growth index 7. We would like to stress 
that the 7-parameterisation of modified gravity has 
its limitations and other more ge neral parameterisa- 



tions h ave been proposed (se e e.g. ISilvestri fc Trodden 
2OO9I : iBean fc Tanamatithaml bOld: IPaniel fc Linde ' 



20ld: IClifton et al.l I2OIII : iHoiiati. Pogosian fc Zhaol |20~ 
Baker et al.ll201ll l. However, this is going beyond the scope 
of this paper. 

We combine our r esult for fle(zefl) wit h the latest 
results from WMAP7 (|Komatsu et all boill ). where we 
use the WMAP7.2 dataset provided on the NASA web- 
pag(Q. While ACDM predicts 7 « 0.55, alternative the- 
ories of gravity deviate from this value. One example of 
su ch an alternative model is the DG P braneworld model 
of iDvali. Gabadadze fc Porratil (|200(]| '). in which our observ- 
able Universe is considered to be a brane embedded in a 
higher dimensional bulk space-time and the leakage of grav- 
ity force propagating into the bulk can lead to the current ac- 
celerated expansion of the Universe. Because of the missing 
dark energy component, thi s m odel predicts a larger growth 
index of 7 0.69 (Lindcr 20051 '). 

In Figure [TS] we compare measurements of the growth 
of structure ge from different galaxy redshift surveys. For 



the WiggleZ survey we include data points which assume 
a correct fiducial cosmology (solid black triangles) as well 
as data points which account for the Alcock-Paczynski ef- 
fect (empty black triangles). The degeneracy between the 
Alcock-Paczynski effect and the linear redshift-space distor- 
tion signal increases the error by about a factor of two. As 
we showed in section lSTSl the Alcock-Paczynski effect is very 
small in 6dFGS, which therefore yields a direct measurement 
of the redshift-space distortion signal. 

All data points seem to be in good agreement with the 
ACDM model (black solid line), while the DGP model gen- 
erally predicts smaller values of gg. The value of as for the 
two different models has been derive d from the CMB scalar 
amplitude As (jKomatsu et al.ll201 j ). where we use the cor- 
responding Friedmann equation to calculate as. 

The analysis method we apply in this section is sum- 
marised in the following four points: 

(i) We produ ce a Monte Carlo Ma rkov Chain (MCMC) 
with CosmoMC (|Lewis fc Bridlell2002t ) for a ACDM universe 
by fitting the WMAP7 dataset. The CMB depends on dark 
energy through the distance of last scattering and the late- 
time Integrated Sachs- Wolfe (ISW) effect. We avoid the con- 
tributions of the ISW effect by limiting the WMAP7 dataset 
to multipole moments £ > 100. 

(ii) Now we importance-sample the CosmoMC chain by 
randomly choosing a value of 7 in the range 7 1 for 
each chain element. Since the value of as{zcf{) depends on 
7 we have to recalculate this value for each chain element. 
First we derive the growth factor 



-D(acff) 



exp 



da f{a)/a' 



(46) 



where a^s is the scale factor at the effective redshift a^H — 
1/(1 -I- Zes)- In order to derive o's,7(-Zeff) we have to extrap- 
olate from the matter dominated region to the effective red- 
shift, 



-(Zcff) = —fr, ^Oi{zm) 



(47) 



50, well in the 



0(zm) 

where we use (Tg(zhi) from eg. 1451 and zy^i 
matter-dominated regime. 

(iii) We now calculate the growth rate using /^(zcff) — 
r2;;,(Zeff) and construct gB,^,{^as) = /7 (zoh) 0-8,7 (2;off). 

(iv) Finally we compare the model with gg — 0.423 ± 
0.055 from Table [T] and combine the likelihood from this 
comparison with the WMAP7 likelihood. 

The result is shown in Figure 1131 Marginalising over 
the remaining parameters we get 7 = 0.547 ± 0.088 and 
Qm ~ 0.271 ± 0.027, which is in agreement with the predic- 
tion of a ACDM universe (7 « 0.55). Our analysis depends 
only on the growth rate measured in 6dFGS and WMAP7. 
This makes our measurement of 7 independent of systematic 
effects like the Alcock-Paczynski distortion which is a matter 
of concern for galaxy redshift surveys at higher redshift. 



FUTURE LOW-REDSHIFT GALAXY 
SURVEYS: WALLABY AND TAIPAN 



In this section we make predictions for the accuracy of fas 
|http: //lambda. gsfc.nasa.gov/product/map/dr4/llkelihood_get . cla^ asurements from future low-redshift galaxy surveys using 
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WMAP7 

6dFGS + WMAP7 
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Figure 13. Likelihood distribution of 7 — Q,m for our fit to 
qg = 0.423 ± 0.055 from 6dFGS and WMAP7 jKomatsu et al.l 
I2OI . Ttie 6dFGS contours (black) include a normalisation prior 
from WMAP7. Marginalising over the different parameters in the 
MCMC chain gives = 0.271±0.027 and 7 = 0.547±0.088. The 
low redshift of 6dFGS makes this measurement particularly sen- 
sitive to 7 and independent of systematic effects like the Alcock- 
Paczynski distortion. 



a Fish er matrix analysis based on IWhite. Song fc Percivall 
(|2008l ) 

The Wide-field ASKAP L-band Legacy All-sky Blind 
surveY (WALLABYfl is an HI survey planned for the Aus- 
tralian SKA Pathfinder telescope (ASKAP), currently un- 
der construction at the Murchison Radio-astronomy Obser- 
vatory (MRO) in Western Australia. The survey will cover 
75% of the sky and a proposal exists to fill up the remaining 
25% using the Westerbork Radio Telescope. I n this analysis 
we fol low the survey parameters employed bv lBeutler et al.l 
(|201ll ) (see also Duffy et al., in preparation), where for WAL- 
LABY we assume a 47r survey containing 600 000 galaxies 
at a mean redshift of 2 = 0. 04. The line ar bias of a typi- 
cal WALLABY galaxy is 0.7 (iBasilakos et al. 2007) and the 
volume of the survey is 0.12/i~"^Gpc^. 

The TAIPAN survejQ proposed for the UK Schmidt 
telescope at Siding Spring Observatory in New South Wales 
will cover a similar sky area as 6dFGS but will extend to a 
larger redshift such that 2 = 0.08. For our TAIPAN forecast 
we assumed the same sky-coverage as 6dFGS (/sky = 0.41), 
a bias of 6 = 1.4, a total of 400 000 galaxies and a volume of 
0.23ft"^ Gpc^ 

First we test the Fisher matrix prediction for /erg in the 
case of 6dFGS. We assume a survey volume of 0.08/i^'* Gpc'^, 
with 81 971 galaxies. Using fcmax ~ O.l/i Mpc~^ we forecast a 
measurement of /erg of 23%, while using fc^ax = 0.2/i Mpc~^ 
produces a 8.3% error. The actual error in /erg we found in 



^ http:/ /www. atnf.csiro.au/roscarch/WALLABY 

* TAIPAN: Transforming Astronomical Imaging surveys through 

Polychromatic Analysis of Nebulae 



this paper is 13%, somewhere between these two values. For 
WALLABY and TAIPAN we will report constraints for both 
fcmax = 0.1/iMpc"^ and fc^ax = 0.2/1 Mpc"\ 

With the specifications given above and using fcmax = 
0.1 (0.2)/i Mpc~\ the Fisher matrix forecast for WALLABY 
is a measurement of /ag with 10.5 (3.9)% error. We in- 
cluded the WALLABY forecast with a 4% error-bar in Fig- 
ure [121 The model TAIPAN survey produces forecast errors 
of 13.2 (4.9)%, improving the results from 6dFGS by almost 
a factor of two. Although TAIPAN maps a larger volume of 
the Universe compared to WALLABY, it does not produce a 
better measurement of /erg. WALLABY has a smaller galaxy 
bias, which increases the redshift-space distortion signal by 
a factor of two compared to TAIPAN. This will also be very 
useful for breaking the degeneracy between bias and erg using 
the technique of section [6l2] 

The WALLABY survey will target galaxies rich in HI 
gas. Such galaxies will mostly populate under-dense regions 
of the Universe, because in groups and clusters galaxies are 
stripped of their gas by interactions with other galaxies and 
the intra-group and intra-cluster medium. This is the reason 
that Hl-selected galaxies possess a low bias (~ 0.7). How- 
ever, this fact also implies that WALLABY galaxies sample 
the density field in a manner that avoids high-density re- 
gions. These high-density regions are an important source 
of non-linear redshift-space distortions ( "finger-of-God" ef- 
fect). We can hence suppose that non-linear effects will be 
smaller in amplitude in an HI survey compared to hiehly- 
biased surveys such as 6dFGS (see e.g. ISimpson et al.ll201ll 
or Figure 4 in lReid fc Whitell201ll ). This should allow the 
inclusion of much smaller scales in the analysis, producing 
more accurate measurements. A more detailed analysis us- 
ing WALLABY mock ca t alogue s is in preparation. 

[McDonald fc SeliakI (|2009l ) have suggested that mul- 
tiple tracers within the same cosmic volume can be used 
to reduce the sampling variance and improve cosmological 
parameter constraints (see also ISeliakI l2009l : ISlosaij l2009l : 
iBernstein fc CailbOllh . Using the ratio of the perturbation 
amplitudes of two surveys with different bias factors gives 



hi + f^? ^ ab2 + /m' 
b2 + /m^ b2 + f^i^ 



(48) 



where a = 61/62. The angular dependence of this expression 
allows one to extract / wi thout any dependence on the den- 
sity field (see Figure 1 in IBernstein fc Gail 120111 ). The den- 
sity field is the source of the sampling variance error, since 
it will change, depending on the patch of the sky which is 
observed. Using the ratio of two tracers, the precision with 
which the growth rate / can be determined is (in princi- 
ple) only limited by the shot noise, and not by the sampling 
variance. 

The three surveys discussed above, 6dFGS, WALLABY 
and TAIPAN, have a large overlapping volume which al- 
lows the potential application of this method. The technique 
works best for densely-sampled surveys with very different 
bias factors, hi and 62. While TAIPAN and 6dFGS have very 
similar bias, the bias of WALLABY will be much smaller. 

We assume an overlap volume of 0.41 x 0.12fe~^ Gpc'^ — 
0.049/i~^ Gpc^, where we multiply the sky coverage of 
6dFGS and TAIPAN with the effective volume of WAL- 
LABY. For the different surveys we use the parameters as 
stated above. We forecast Fisher matrix constraints on /erg 
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of 10.3 (5)% using fc^ax = 0.1 (0.2)/iMpc-^ Because the 
overlap volume is only ~ 1/3 of the WALLABY volume, 
this result does not improve the measurement arising from 
WALLABY alone, especially considering that WALLABY 
itself may be able to include modes up to large fcmax in the 
fitting process. 

Our results show that future surveys such as WAL- 
LABY and TAIPAN will provide an accurate measure- 
ment of /erg at low redshift, and will be able to com- 
plement future high-redshift surveys such as BOSS, which 
will have a similar accuracy for se veral data points over 



the higher redshift range 0.2 - 0.6 iS ong & Percivall [20091 : 

' 31: 



IWhite. Song fc PercivallboOSl : iReid fc Whitc..201lh . 



• We used a Fisher matrix analysis to forecast the con- 
straints on fas that would be obtained from two future 
low-redshift galaxy surveys, WALLABY and TAIPAN. We 
found that WALLABY will be able to measure /ag to a 
forecast accuracy of 10.5% for fcmax = O.l/iMpc"^ and 3.9% 
for fcmax = 0.2/1 Mpc"\ A combination of 6dFGS, TAIPAN 
an d WALLABY, using t he mu ltiple-tracer method proposed 
bv lMcDonald fc SeliakI (|200^ ), will be able to constrain /(Tg 
to 5 — 10.3%. These measurements would complement fu- 
ture large-volume surveys such as BOSS, which will mea- 
sure the growth rate at much higher redshift [z > 0.2), and 
contribute to future precision tests of General Relativity on 
cosmic scales. 



9 CONCLUSION 

In this paper we have measured the 2D correlation func- 
tion of the 6dF Galaxy Survey. We derived a covariance 
matrix using jack-knife resampling as well as log-normal re- 
alisations and showed that both techniques give comparable 
results. We have modelled the 2D correlation function with 
a simple streaming m ode l and a more advanced approach 
suggested bvIScoccim arrd (12004) combined with the N-bod y 
calibrated results from Jennings. Baugh fc PascolU (|2011al ). 
We formulated these models in real-space including wide- 
angle cor rections. For t he fin al results on /ag we chose the 
model by [Scoccimarr 

3 (|2004 ). although we found that both 
models gave consistent results at sufficiently large scales. 

We analysed the measurement in two different ways. 
First we fitted for the two parameters g0{zcs) — 
f {zcft)as{zctf) and gbizcs) = 6cr8(2:off), where these con- 
straints depend only on the GdFGS data. Our second anal- 
ysis method assumes a growth index from standard gravity 
(7 « 0.55) and fits for ag, b, Q,n, Hq and ap, where we com- 
bine the GdFGS measurement with a prior in the Hubble 
constant. All parameter measurements are summarised in 
Table [T] We can summarise the results as follows: 

• Our first analysis method found gg{zcfi) = 
/(zcff)o"8(2off) = 0.423 ± 0.055 and gb{zcfi) = basiz^a) = 
1.134 ± 0.073, at an effective redshift of z^a = 0.067. 
The 6dFGS measurement of ge, unlike high-redshift mea- 
surements, does not depend on assumptions about the 
expansion history of the Universe and the Alcock-Paczynski 
distortion. 

• In our second analysis method we used the angle de- 
pendence of redshift-space distortions in the 2D correlation 
function to break the degeneracy between the galaxy bias 
6 and the normalisation of the matter clustering statistic 
(Tg, assuming standard gravity. We found erg = 0.76 ± 0.11, 

= 0.250 ±0.022, b = 1.48 ±0.27 and = 174±73 km/s. 
This result uses a prior on Ho from lRiess et al.l (|201ll ). 

• Gombi ning our measurem ent of gg{zcB) with 
WMAP7 (|Komatsu et all I2OIII ') allows us to measure 
the growth index 7, directly testing General Relativity. We 
found 7 = 0.547 ± 0.088 and fi™ = 0.271 ± 0.027, in agree- 
ment with the predictions of General Relativity (7 « 0.55). 
The 6dFGS measurement of this parameter is independent 
of possible degeneracies of 7 with other parameters which 
affect the correlation function at high redshift, such as the 
dark energy equation of state parameter w. 
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APPENDIX A: PARTIAL ANALYTICAL 
SOLUTION FOR THE CORRELATION 
FUNCTION MOMENTS INTEGRAL 

The double integrals in eqs. [35] and [38] are difficult to solve 
numerically. Here we show an analytical solution for the in- 
tegral over ^ which allows for a faster numerical solution of 
the full integral. First we re- write eq. [38] as 



.-1 (27r)2 
X cos{kry)Dxy{k)'Pi{^Ji) 



(2^)' 

where i is the complex number. More generally the integral 
over ^ can be written as 



(/i) = / d^e 
. d 



ikrfj, — {kiJ,au ) 



d{kr) 
d 



dfi e 



(Al) 



T./2 



ikr ^ — {kjia-u ) 



The integral on the right can be solved analytically. If we 
set a = kr and b — {ka^)^ we obtain 



^ —ia — b 

Vb 



Daw 



a — 2ib 
2Vb 



2iaj^ I a + 2ib 
e Daw 



2Vb 



(A2) 



Taking the n-th derivatives of the real part of the term above 
gives Fn{^). The Dawson integral Daw(a;) can be calculated 
using the imaginary error function erfi(2;): 



Daw(x-) = e 



e^ dy 



erfi(a;). 



(A3) 



We can than construct eq. [38] for the different correlation 
function moments: 



€,T=2,xy{''') 
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iT=A,xy{'r) 
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2Vx,(fc)-[3F2(/i)-Fo(M)] 



(27r 

°° k'^dk 



2V.!/(fc)2 [5J'3(/i)-3Fi(/i)] 







(^v.!/(fc) g [35F4(/i) - 30F2(/i) + SFoif^)] 



(A4) 

and so on. 
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